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Abstract 

Considering the density wave of scalar and pseudoscalar condensates, we study the response of quark matter to a weak external 
magnetic field. In an external magnetic field, the energy spectrum of the lowest Landau level becomes asymmetric about zero, which 
is closely related to chiral anomaly, and gives rise to the spontaneous magnetization. This mechanism may be one of candidates for 
the origin of the strong magnetic field in pulsars and/or magnetars. 
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Recently, the existence of the inhomogeneous chiral phase in 
the QCD phase diagram has been discussed by the analysis of 
the effective models such as Nambu-Jona-Lasinio (NJL) model 
Qii or the Schwinger-Dyson approach lil]. In this phase, 
the quark condensates spatially modulate and it is very similar 
to the FFLO state in superconductor Hi or spin/charge den¬ 
sity wave llTllst]. Here, we consider “dual chiral density wave 
(DCDW)” yj] among many kinds of form of the condensates: 
the quark condensates then take the form, 

A(r) s {ij/if/) + i{ij/iy^Titf/) - Ae'^^, (1) 


within the two-flavor QCD. This configuration is also obtained 
by embedding one of the Hartree-Fock solutions in the NJL 2 
model, so-called chiral spiral HUS]- Since the DCDW phase 
has been expected to appear in the moderate density region li, 
it may be plausible that this phase is realized in neutron stars. 

The effect of the magnetic field has been first discussed by 
Frolov et al. for the DCDW phase fUll . They have found that 
the spatial direction of the wavevector q is favored to be paral¬ 
lel to the magnetic field, and the domain of the DCDW phase 
is much extended in the QCD phase diagram. In ref. iH these 
features arises from some topological effect through spectral 
asymmetry of the quark energy; quarks exhibit an interesting 
feature in the presence of the magnetic field and the energy 
spectrum becomes asymmetric about zero. There also appear 
new terms in the generalized Ginzburg-Landau expansion due 
to spectral asymmetry, which signals the novel Lifshitz point in 
the QCD phase diagram. Thus, they emphasized the peculiar 
role of the phase degree of freedom of A(r). 

Here we further inquire this issue. We study magnetic prop¬ 
erties of the DCDW phase to reveal another aspect, spontaneous 
magnetization in the DCDW phase, which suggests a micro¬ 
scopic origin of the strong magnetic field in compact stars. 

The origin of the strong magnetic field in compact stars has 
been one of the long-standing problems. In particular, mag¬ 
netars have the huge magnetic field ~ lO'^G on the surface 
liSllldl]- As a candidate of the origin, amplification of the mag¬ 


netic field by the dynamo mechanism, magnetorotational insta¬ 
bility or the hypothesis of the fossil magnetic field has been 
proposed so far from the macroscopic point of view. Although 
numerical simulations have been actively performed, no def¬ 
inite conclusions have been obtained. From the microscopic 
point of view, it has been proposed that the spontaneous mag¬ 
netization emerges by spin alignment of quarks on the analogy 
of the electron gas 11.511 . However, this phase should be devel¬ 
oped in the low density region. As another mechanism, it has 
been proposed that axial anomaly acting on the parallel layer of 
the pion domain wall produces magnetization in nuclear matter 


mil 


We use the two-flavor NJL model in the mean field approxi¬ 
mation. It is sufficient to consider the each flavor case because 
Lagrangian is diagonalized about the flavors. Then, the Dirac 
Hamiltonian takes the form. 


H — -ia ■ D — 2Gy 




( 2 ) 


with the covariant derivative, D = V -H iejA. Taking the ex¬ 
ternal magnetic field B along the z axis, the energy spectrum 
constitutes the Landau levels in, 


^ -H q/lj + 2\efB\n, {n + 0) (3) 

= e y]pz+ + ?/2, (« = 0) (4) 

with /"=+!, which denotes the spin polarization. For m < ql2, 
the sign of e does not always correspond to the particle or anti¬ 
particle state because the lower branch of the lowest Landau 
level (LLL) in - 0) is not always negative. In the higher Lan¬ 
dau levels (hLLs) (n + 0), there are four energy branches. On 
the other hand, LLL has only two energy branches and becomes 
asymmetric about zero. As a result, the thermodynamic poten- 
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where 


tial takes the form in the two-flavor case. 


Q(//, T, B-m,q)^ — + Nc y (5) 


where, 

Q.f = 


+ J]in[ajl + (Ep,e-Mf]i ( 6 ) 


with the Matsubara frequency, lOk - (2k -i- l)7rr. To investigate 
the response of quark matter to the weak magnetic field B, the 
thermodynamic potential is expanded about B, 


T,B ',m,q) — 

Q®(//, T ■,m,q) + eB T ■, m, q) + ■ ■ ■ , (7) 


where e denotes the elementary charge. Since the vacuum part 
in has divergence, it must be regularized by ,e.g., the proper 

time regularization (PTR) llj]. LLL contributes only to Q*'' be¬ 
cause the energy spectrum does not depend on B and the B de¬ 
pendence only emerges through the Landau degeneracy factor. 
On the other hand, hLLs contribute to the all order terms of B. 

The magnetization can be deduced from the thermodynamic 
potential as. 


M = - 


T,B) 

dB 


( 8 ) 


where Qmin represents the minimized thermodynamic potential 
about the order parameters and only depends on /r, T and B. 
Taking the limit, B ^ 0, we find the spontaneous magnetization 
in the form. 


Mo = T\m — q — 


(9) 


where and q^^^ represent the minimal values for In the 
following, we will figure out the peculiar role of LLL and show 
that it leads to spontaneous magnetization. 

For the evaluation of we must carefully treat the effect 
of chiral anomaly. According to Refs. USUI], spectral asym¬ 
metry generally gives rise to anomalous particle number. 




lim ^ y sign(/l^)|T,t| ", 

9^+0 Z Z—^ 


( 10 ) 


where is the eigenvalue of the arbitrary Dirac Hamilto¬ 
nian. Spectral asymmetry is ill-defined as it is and needs a 
proper regularization without violating the gauge invariance. In 
the DCDW phase, LLL exhibits spectral asymmetry to induce 
anomalous particle number proportional to B OITII . Then, the 
LLL contribution in can be decomposed into three terms. 


^(l).LLL ^ ^(1) LLL ^ ^(1),LLL ^ ^(D.LLL 


( 11 ) 




+ 



^a).LLL ^ 


NcT 

In 


f 




( 12 ) 


(13) 

(14) 


with Wf = e + nP- + qjl. The density dependent term 

Q^i)’Lll includes the anomalous contribution, caused by 

spectral asymmetry. The //-invariant, qn, renders 


qn = 


lim fpy 

.s->+o J 2 ti Zj 


|Wf| ■'sign(Wf) 



(m > ql2) 
(m < ql2) 


(15) 


When m > ql2, this quantity agrees with the contribution 
of the chiral anomaly represented by the Wess-Zumino-Witten 
(WZW) term tH. The WZW term does not depend on m but it 
vanishes in the limit, m —> 0. The contribution of hLLs to 
should be carefully evaluated by expanding the thermodynamic 
potential with respect to B after the summation over n. Then, 
the hLLs contribution in can be similarly decomposed into 
three terms. 


Q(l),hLL ^ Qa),hLL QillhLh ^ Qa).hLL^ 


(16) 


which does not include the anomalous contribution since hLLs 
have no spectral asymmetry. We find that the three terms are 
the even function of q, and Thus = 

q(1),lll Q(i),hLL jjQj (jiverge without any regularization 
and renders 




1^ 

An 

An 


qn 


ffzz 


t(/LI - TOjP) 0 (t(jJ^) 0 (b - TUe) 


- 07, 

e r=±l 


The first term can be interpreted as the contribution of anomaly 
and the second and third terms as the contribution of positive- 
energy valence quarks. Note that the even function of q in 
Eq. (fTTl) is completely canceled by the corresponding one in 
Eq. (fTST l to make the odd function of q. It vanishes in the 
limit, m —» 0, which behavior may be physically reasonable 
because there should be no condensate and q should be redun¬ 
dant in this limit. Consequently we can see that emerges 
only when m ^ 0,q ^ 0. In other words, quark matter has the 
spontaneous magnetization only in the DCDW phase. 

Note that the spontaneous magnetization takes a special form 
for some peculiar values of the parameters, m, q, as well as p, T. 
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{a)T = 0 (b) r ^ 30MeV 

Figure 1: The chemical potential dependence of the order parameters and spontaneous magnetization at B = 0 


At ju < m - qjl and T — 0, where is no valence quark, the mag¬ 
netization only comes from the WZW term, Mq - argued 
in Ref. ifl^ . Unfortunately, such a situation is not realized 
in the present calculation because there always exist valence 
quarks in the DCDW phase, which is given rise to by the nest¬ 
ing effect of the Fermi surface Oil ItL l20l] . On the other hand, 
when m becomes small compared to p or F but q still is not 
small. Mo is evaluated to be 


Mo = 


Nc r.. /I 

-<■&? Z 5 ■"' 

cr=±l ' 

0{m^^ , 


.ql2-o-fi\ 2 
m 


2nT 


(18) 


where if/ is the digamma function. Then, the leading order about 
m is by the symmetry, m —> -m. 

The magnetization is numerically evaluated at T = 0 and 
T ^ 30MeV, for example. In this calculation, we choose 
A = 660MeV as the cutoff in PTR and GA^ = 6.35, which 
reproduce pion decay constant f„ - 93MeV and constituent 
quark mass ^ 330MeV in the vacuum. FIG.[T]shows the chem¬ 
ical potential dependence of the order parameters and the spon¬ 
taneous magnetization Mq. There are three phases, the DCDW 
phase (m 0, q + 0), the homogeneously chiral-broken phase 
(m ^ 0, q - 0) and the chiral-restored phase {m = 0). We 
can see that Mq becomes nonzero only in the DCDW phase 
and has discontinuity at the transition point from the homoge¬ 
neously chiral-broken phase. At this point, the order parameters 
also discontinuous, which implies the phase transition is of the 
first order. On the other hand, near another transition point to 
the chiral-restored phase, the spontaneous magnetization is pro¬ 
portional to Furthermore, as temperature increases, the 

spontaneous magnetization decreases and the region of magne¬ 
tized phase gets narrower. 

Considering a sphere of the uniform magnetization Mq, the 
magnetic field, B - yMo is produced on the surface. For the 
sphere of quark matter with a constant density in the DCDW 
phase at zero temperature, the magnetic field is estimated B ~ 
lO'^G on the surface, which may be comparable with the obser¬ 


vation of magnetars. However, this estimate might be too rough 
because density is not constant in neutron stars. Finite shell 
structure of the DCDW phase may be another possible realiza¬ 
tion. Considering that the DCDW phase gets narrow at high 
temperature, magnetic field may become much lower than this 
estimate in hot neutron stars. Actually the spontaneous magne¬ 
tization vanishes with the disappearance of the DCDW phase at 
T ~ lOOMeV. At r ^ 30MeV, which may be a relevant tem¬ 
perature in newly born neutron stars, the magnetic field remains 
to be a similar order of magnitude. 

In Ref. Ill], the expectation value of the magnetic moment, 
(^(r)cr'^i/f(r)), is evaluated, which behaves like spin density 
wave and vanishes after the spatial average. This seems to con¬ 
tradict with the present results. However, we can see that a 
careful application of the Gordon identity leads to the present 
results. The external magnetic field minimally couples with the 
quark field through the covariant derivative, J d‘^x in 

the Lagrangian, where Q is the electric charge matrix in the fla¬ 
vor space. According to Ref. 12 111 , this term can be decomposed 
into the form. 


/ 




- f 

2m J 


cfx 


xQo-^^X +xQ'^ixD2X 


+ xQcr^^i^Tiqxx , 


(19) 


where the Landau gauge is taken, A = (0, Bx, 0), and x rep¬ 
resents the quark held after the Weinberg transformation, x - 
gir Here we have used the modified Gordon decompo¬ 

sition by using the Dirac equation in the presence of DCDW, 
instead of the usual treatment by the free Dirac equation. Con¬ 
sequently, we have the expectation value of the magnetic mo¬ 
ment, Ct(r)cr^^;^f(r)), instead of the above one in terms of ij/, 
which never vanishes after the spatial average. There appear 
two contributions in ([19]) besides the magnetic moment: the 
second term can be interpreted as the angular momentum and 
the third term proportional to q comes from the operator inher¬ 
ent in the DCDW phase. Thus, the magnetization discussed 
here should be regarded as the statistical average of these oper- 
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ators. 

We have shown that quark matter has spontaneous magne¬ 
tization in the DCDW phase where the response to the exter¬ 
nal magnetic field changes from the one in the homogeneously 
chiral-broken phase or the chiral-restored phase, that is, 
emerges in the thermodynamic potential. In Ref. in, they 
have simply discussed the spontaneous magnetization by using 
the WZW term in the effective Lagrangian written in terms of 
only mesons. However, we have seen that has not only 
the anomalous contribution closely related to the WZW term 
but also the contribution of valence quarks. In this sense, the 
fermion degrees of freedom and their dynamics are indispens¬ 
able. Thus we can say that DCDW provides a realistic ground 
where the spontaneous magnetization is realized. 

The spontaneous magnetization discussed here has an inter¬ 
esting feature different from the usual spin alignment iflfill : the 
magnetization in the DCDW phase is caused by the different op¬ 
erators from the naive magnetic moment. Further discussions 
about the peculiar magnetic properties such as magnetic sus¬ 
ceptibility or the Nambu-Goldstone mode (spin wave) will be 
presented in another paper. 

We have emphasized that the spectral asymmetry is im¬ 
portant for the mechanism of the spontaneous magnetization. 
Since the complex order parameter^r) is necessary for the 
energy spectrum to be asymmetric lll2|] . it is conceivable that 
there is no spontaneous magnetization in the “real kink crystal” 
phase with the real order parameter A(r). On the other 
hand, the spontaneous magnetization should emerge in the sim¬ 
ilar way in the phase with hybrid chiral condensate 11221] . 

It has been discussed that the NIL model, in the mean field 
approximation, exhibits magnetic catalysis 123, mm, while 
the recent lattice calculation shows inverse magnetic catalysis 
1 2^ 271. To explain this phenomenon an idea of the effective 
four-Fermi coupling has been proposed within the NIL model 
1 2a 1291], where the coupling strength depends on the external 
magnetic field by considering the coupling of the magnetic field 
with the quark and gluon loops constituting the four-Fermi cou¬ 
pling. Even if this is the case, however, the response to the tiny 
external magnetic field may be little changed and our findings 
should not be affected. 

In order to explain the strong magnetic field in neutron stars 
by the present mechanism, it is necessary to study the appear¬ 
ance of the DCDW phase in more realistic conditions such as 
the charge neutrality or the chemical equilibrium. It is also an 
open question that the DCDW phase survives when quarks have 
the finite current mass M- 
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